Probing the homogeneous spectral function of a strongly interacting superfluid atomic Fermi gas in 
a trap using phase separation and momentum resolved rf spectroscopy 
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It is of central importance to probe the local spectral function yl(k,aj) of a strongly interacting Fermi gas 
in a trap. Momentum resolved rf spectroscopy has been demonstrated to be able to probe the trap averaged 
A(k,cj). However, the usefulness of this technique was limited by the trap inhomogeneity. Independent of a 
specific theory, here we propose that by studying the momentum resolved rf spectra of the minority fermions of 
a phase separated, population imbalanced Fermi gas at low temperature, one can effectively extract j4(k, uj) of 
a homogeneous superfluid Fermi gas (at the trap center). In support, we present calculated spectral functions 
and spectral intensity maps for various cases from BCS through BEC regimes using different theories. 

PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.-z 
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I. INTRODUCTION 

Ultracold atomic Fermi gases has emerged as a rapidly de- 
veloping field, bridging condensed matter and atomic physics. 
Owing to their easy tunabilities such as the effective two-body 
interaction strength, population imbalance, and dimensional- 
ity, atomic Fermi gases can be viewed as a quantum simula- 
tor of many important condensed matter systems, such as the 
Hubbard model and high Tc superconductors |lll|31- 

Of central importance in a many-body system is the single 
particle spectral function A(k, ui), especially in the strongly 
interacting regime (where different theories often do not 
agree); it can be used to calculate essentially all other physi- 
cal quantities, and to test various theories. For example, one 
important question, which has been long standing for high Tc 
superconductors, is whether a pseudogap exists and how it 
evolves in the superfluid state. Because measurements often 
involve integration of the entire momentum space and the en- 
tire trap, the intrinsic spatial inhomogeneity introduced by the 
trapping potential constitutes a severe problem. It has been a 
long-sought goal to extract the local or homogeneous proper- 
ties, especially A{k, lu), from measurements in the trap. 

In this paper, we propose that one can use momentum re- 
solved rf spectroscopy in conjunction with a high popula- 
tion imbalance to measure the homogeneous spectral function 
A(k, oj), for the density and chemical potential at the trap cen- 
ter. 

Since the first experimental realization of BCS-BEC 
crossover, the most direct experimental probe of the single 
particle properties is arguably the rf spectroscopy. The first 
generation of rf measurements involves integration over the 
entire trap and the whole momentum space which led to 
a two-peak structure of the rf spectra at an intermediate tem- 
perature T below Tc. Due to the integrations, controversies 
arose regarding the origin of this two-peak structure. Tomo- 
graphic rf spectroscopy 0] involves integration over the entire 
momentum space, and therefore cannot be used to probe the 
spectral function. The recent Ho-Zhou scheme |5] allows one 
to calculate the density and fermionic chemical potential ji as 
a function of the radius r. However, it relies on the assumption 



of a strict harmonic trapping potential f6'l, often not satisfied. 
More importantly, it does not provide any microscopic infor- 
mation such as the spectral function. 

The controversies regarding the origin of the aforemen- 
tioned double peak structure have largely been cleared by the 
recent experiments |7] on using momentum resolved rf 
spectroscopy. With momentum resolution, the rf spectroscopy 
essentially measures the centrally important A(k, w), which is 
of central importance in a many-body system. However, this 
is true only for a homogeneous Fermi gas. The inherent spa- 
tial inhomogeneity severely limits the quantitative resolution 
of the extracted spectral function in the experiment of Ref. 
which involves integration over the entire trap. 

Here we address this inhomogeneity issue and show that it 
can be largely avoided utilizing phase separation at high pop- 
ulation imbalances. In what follows we will use two differ- 
ent theories to do computations. Despite some differences in 
quantitative details, we emphasize that the validity of our pro- 
posal is model independent and can be applied in other theo- 
ries as well. Most importantly, it is a proposal for experiment. 

As mentioned elsewhere |i8J, the use of '"^K rather than ^Li 
is crucial in extracting A(k, w), since there are no complica- 
tions from final state interactions l'9'l near the usual Feshbach 
resonance around 202 G. It is for this reason that the rf spectral 
intensity for is simply proportional to the spectral func- 
tion A(k, cj). For ^Li, extra efforts are needed in order to 
extract ^(k, a;) from the rf spectral measurements. 

The spectral function v4(k, cj) can be used to uniquely con- 
struct the Green's function and the single particle self energy, 
for which different theories often give different results. There- 
fore, the more accurate the measurement of ^(k, lo), the eas- 
ier it is to test different theories. For example, different BCS- 
BEC crossover approaches ifTol [Till show different dispersive 
behavior of yl(k, w). The theoretical scheme used in Ref. 
exhibits a clear downward dispersion in the spectral intensity 
map near Tc, manifesting an existing pseudogap before the 
onset of superfluidity at lower T — Tc, different from those 
[12, 13] that follow the approach of Nozieres and Schmitt- 
Rink (NSR) |14]. A quantitatively accurate measurement of 
A{k,ijj) should serve to unambiguously test these different 
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schools of theories. 



II. THEORETICAL FORMALISM 

In a typical rf spectroscopy measurement, pairing takes 
place between two low-lying hyperfine states, which we re- 
fer to as levels 1 and 2. An rf field of frequency v is used 
to excite the atoms in hyperfine level 2 to another hyperfine 
state, which is unoccupied initially and referred to as level 3. 
It has been shown previously [8J that for ''^^K we have for the 
momentum-resolved RP current 
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where A(k, lS) — —2 Im G(k, a; + iO), ](x) the Fermi distri- 
bution function, with the transition matrix element set to unity. 
Here G is the Green's function of level 2 atoms, which we take 
to be the spin down or minority species. Similar to angle- 
resolved photoemission spectroscopy in a usual condensed 
matter system, the rf current measures v4(k, lS) directly. Note 
that here bv is the rf detuning. The frequency uj — (^\^ — 8v 
corresponds to the energy of level 2 atoms measured with re- 
spect to their Fermi level, where = fc^/Sm — [i, [iv& the 
chemical potential of atoms in level 2. We take h ^ 1. It 
is obvious that when level 2 atoms are free, we have — 
and u! — the former gives the sharp peak at zero detun- 
ing in previous rf spectra in Ref. '3! whereas the latter gives 
the free atom dispersion in the spectral intensity map in the 
(jj — k plane in Ref.j?! As has been used in the experiment |7], 
the angle-integrated "occupied spectral intensity" is given by 
IP^°*°{k,uj) = {k'^/2'K^)A{k,oj)f{uj). 

The central issue here is to calculate the spectral function 
or equivalently the Green's function G{K) for the minority 
(level 2) atoms. Here we present calculations using two dif- 
ferent approaches. In the first, pairing fluctuation, approach, 
in the presence of population imbalance, detailed calculations 
of the superfluid phase diagram, the Green's function and den- 
sity profiles in each phase can be found in Refs. ^andHS In 
the case of phase separation in a trap, as shown in Fig. 2 of 
Ref. OH there is a BCS superfluid core (without population 
imbalance), surrounded by the majority atoms, as confirmed 
experimentally ITtI [Tsll . To a first order approximation, the 
density profile of the minority atoms can be obtained by trun- 
cating the density profile of an unimbalanced Fermi gas at the 
(sharp) phase separation boundary. With and without popu- 
lation imbalance, the normal state self energy Spg follows a 
rather simple BCS-like form [8, 19J. In the superfluid state, 
the self energy I^{K) contains two terms, associated with the 
condensed (Use) and noncondensed fermion pairs (Epg), re- 
spectively; I](k, oj) = I]pg(k, w) + I]sc(k, w), where 
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Here Apg, A^^ and A — \J + A^^ are the pseudgogap, 
superfluid gap and total excitation gap, respectively. The 
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Figure 1. (Color online) Radial profile of the excitation gap A/E% 
in the superfluid phase of a trapped Fermi gas without population 
imbalance for l/k%a = —0.5, 0, 0.5 and 1, corresponding to BCS 
through BEC regimes, calculated at T = 0.5rc, using self energy 
Eq. The gap is nearly flat for r < O.SRtf- 



broadening 7 7^ and "incoherent" background contribu- 
tion So can be determined by fitting the experimentally mea- 
sured rf spectra. The precise values of 7 and 'Sq, and their 
T-dependencies are not particularly important for the present 
purposes. Although the pseudogap contribution becomes less 
pronounced when the temperature is low enough to exhibit 
phase separation, it still makes a significant difference in com- 
parison with a strict BCS theory. 

In the second, mean-fleld, approach, as a test case, we use 
the same density profiles as obtained from the first approach, 
but take a simple broadened BCS self energy. 
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Both approac hes give th e same form of quasiparticle disper- 
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sion = + A^ albeit with different meanings of A. 



Obviously, they become equivalent above T^, when we take 
7 = Eq. In contrast, they differ dramatically below Tc ifioll 
because one contains a pseudogap while the other does not. 
In both approaches, 7 and Sq scale with local Ep{r). As in 
Ref. 8, here we will take Sq as T independent, 7 linear in 
T/Tc, where Tc is calculated using the first approach. The 
rf spectra are finally convoluted with a Gaussian broadening 
function with a standard deviation a, reflecting the instrumen- 
tal resolution in experiment caused by e.g. finite energy and 
momentum resolution of both the rf pulse and the time of 
flight imaging technique. 



III. RESULTS AND DISCUSSIONS 

We first study the radial profile of the excitation gap in 
the superfluid phase. While throughout this paper we study 
each case in a harmonic trap with a local density approxima- 
tion using both approaches, here we show in Fig.[T]the result 
from the pairing fluctuation approach Eq. (|2]). It suffices to 
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Figure 2. (Color online) Occupied spectral intensity maps 
I'''^°*°{k,u}) for the minority atoms of a population imbalanced, 
phase separated unitary Fermi gas in a harmonic trap for (a) R^ = 
O.IRtf, (b) 0.2Rtf, (c) O.SRtf, and (d) OARtf, respectively, 
calculated at T = 0.13T]?, using self energy Eq. The (magenta) 
dashed curves represent the local (or homogeneous) dispersion at the 
trap center, and the (white) solid curves are the quasiparticle disper- 
sion given by the loci of the peak location of the EDC. The broaden- 
ing parameters Eo — 7 (7c) = 0.1 Ep at the trap center. Here the 
spectral intensity increases from (dark blue) continuously through 
the maximum (dark red). The quantitative scale for, e.g., panel (c) 
may be seen from Fig. |4] For the instrumental Gaussian broaden- 
ing, we take a = 0.2Ep (as extracted |8] from the experiment in 
Ref.0). Up to Ri = OARtf, the two sets of curves are essentially 
indistinguishable. 



calculate for cases without population imbalance, since upon 
phase separation the superfluid core becomes unimbalanced 
so that the gap profile of minority atoms can be obtained by 
cutting off the curves in Fig. [1] at different radii for differ- 
ent population imbalances p = {N^ — Ni) / {Nf + Ni). For 
simplicity, here the data were all calculated at T = O.STc, 
below which phase separation takes place at unitarity ifTsll . 
The curves, as labeled, correspond to l/k^a = —0.5 (BCS 
case), (unitary), 0.5 (pseudogap) and 1.0 (BEC) cases, re- 
spectively, where a the inter-fermion s-wave scattering length, 
E^p = keT^ = T?{k%fl2m and k% are the global Fermi en- 
ergy and Fermi momentum in the noninteracting limit for the 
majority atoms. Figure[T]reveals that for r < 0.3i?TF, the gap 
is nearly flat from BCS through BEC regimes, where Rtf is 
the Thomas-Fermi radius. Since the pseudogap increases with 
pairing strength and it is already comparable with the zero 
T gap at unitarity 13, 01. the total gap profile A(r) remains 
nearly flat for small r even at higher T « T^, at unitarity and 
in the BEC regime S. 

With this knowledge, now we study the occupied spec- 
tral intensity maps of the minority (spin down) atoms in the 
(w + \x[r)) - k plane for different degrees of population im- 
balance p, as characterized by the radius i?^ of the minor- 
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Figure 3. (Color online) Occupied spectral intensity maps 
P^°^°(k,iJ) for the minority atoms of a population imbalanced, 
phase separated Fermi gas with R^ — 0.3 Rtf for (a) l/k%a = 
—0.5 (BCS case), (b) (unitary), (c) 0.5 (pseudogap), and (d) 
1.0 (BEC case), calculated at T = 0.09r|,, 0.13r|,, 0.15T^ and 
0.15Tp, respectively, assuming self energy Eq. ([3}. The legends are 
the same as in Fig. |2] The parameters Eo = 7(Tc) at the trap cen- 
ter are (a) 0.1E%, (b)-(c) 0.2E%, and (d) CSSB^, respectively. Tc is 
taken from the pairing fluctuation approach. The instrumental broad- 
ening a — 0.2Ep. For all cases, the two sets of curves are essentially 
indistinguishable. 



ity atomic gas. For illustration purpose, we present the result 
from the pairing fluctuation approach Eq. (|2]i for the unitary 
case with phase separation in Fig. |2] for different R^. Note 
here that, as in Ref. 8, we use ut + /i(r) instead of uj in the 
vertical axis, since the former combination is independent of 
r representing the single particle energy measured from the 
bottom of the band. While Fig. |2a) corresponds to an ex- 
tremely high population imbalance, which may not be readily 
realizable in experiment, Fig.|2c)-(d) are certainly accessible 
experimentally |17]. For example, for ^Li, Ref. lilSil reports 
a phase separation boundary at about O.SRtf at p = 0.54. 
The (magenta) dashed curves represent the local (or homo- 
geneous) dispersion at the trap center, and the (white) solid 
curves are the quasiparticle dispersion given by the loci of the 
peak location of the energy distribution curves (EDC). It is 
evident that up to R^ ~ Q.4Rtf, the two sets of curves are 
essentially indistinguishable. 

Shown in Fig. [3] are the occupied spectral intensity maps 
of the minority atoms of population imbalanced, phase sep- 
arated atomic Fermi gases in a harmonic trap with minority 
radius R^ = Q.3Rtf for different pairing strength from BCS 
through BEC. Here we show the results calculated using the 
mean field self energy Eq. Q. Since the atomic cloud shrinks 
with increasing 1/kpa, these same R^ corresponds to an im- 
balance of about p = 0.9 ^ 0.7. Here the parameters 7 and 
So increase from BCS to BEC, reflecting an increasing exci- 
tation gap. Clearly, for all cases shown, the dispersions ex- 
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Figure 4. (Color online) Comparison of occupied spectral func- 
tion A{k,Lo)j{ijj) between atoms at the trap center (black solid 
and red dashed curves) and trap averaged minority atoms (green 
long dashed and blue dot dashed curves), for a population imbal- 
anced, phase separated unitary Fermi gas with = 0.3-Rtf at 
T = O.STc = O .lSrg. The spectra are calculated for fc = 0.88fc^ f» 
near the Fermi surface, using self energy Eq. (|2j. 
The narrow peaks are the intrinsic spectral line, whereas the broad 
ones are those convoluted with a instrumental broadening function. 
The broadening parameters Eo, 7 and a are the same as in Fig. [T] 
The trap averaged minority spectral function is nearly the same as 
that at the trap center. 



tracted from the trap averaged minority atoms are very close 
to that calculated at the trap center. 

Note that the upward dispersing branch of thermally excited 
quasiparticles is completely suppressed by /(w) and becomes 
invisible here. However, its leftover at large k form a very 
weak broad peak at negative w, as indicated by the light blue 
area above the downward dispersing branch. At higher T, this 
feature ll2lll naturally explains the artificial abrupt jump in the 
quasiparticle dispersions extracted at higher T using a single 
peak fit in Fig. 1 of Refill- 
Finally, shown in Fig.|4]is a direct comparison of the (oc- 
cupied) spectral function between the atoms at the trap center 
(black solid and red dashed curves) and the trap averaged mi- 
nority atoms (green long dashed and blue dot-dashed curves) 
of a highly population imbalanced, phase separated Fermi gas, 
from the pairing fluctuation approach. The curves are for 
k = 0.88fc^ sa y/2mn{r = 0), which is near the Fermi level 
[where A{k,uj) is most sensitive to the variation of A] at the 
trap center The spectral peaks at +Ey^ « A are suppressed 
by /(oj). The narrow peaks are the spectral lines without in- 



strumental broadening, and the broad ones are the same re- 
sult but with instrumental broadening. There are small differ- 
ences between the unconvoluted spectral peaks (black solid) 
and (green long dashed curves) in the peak location and peak 
height (or equivalently peak width). This is primarily due to 
the gradually decreasing broadening parameters caused by de- 
creasing local Ep{r) as r increases, leading to slightly higher 
peak and smaller average gap (given by the peak location). 
We emphasize that these differences are barely discernible and 
are completely washed out by instrumental broadening, as re- 
vealed by the complete overlap between the two broad peaks. 
The complete spectral function A{k, u) can be obtained by 
dividing the unconvoluted lines by /(cj). 

Indeed, a comparison between Fig. |2fc) and Fig. [3jb) re- 
veals that the spectral lines from self energy Eq. ^ have a 
broader background than that from Eq. (|2|. Precise spectral 
function measurement as shown in Fig. |4] is expected to tell 
them apart. 

IV. CONCLUSIONS 

In conclusion, in the presence of phase separation for a 
highly population imbalanced, strongly interacting Fermi gas, 
the trap averaged quasiparticle dispersion of the minority 
atoms is very close to that at the trap center. Therefore, one 
can use the former dispersion to extract effectively not only 
the excitation gap and the fermionic chemical potential but 
also the centrally important spectral function at the trap cen- 
ter This scheme does not depend on our particular BCS-BEC 
crossover theory, nor does it depend on the strict harmonic- 
ity of the trapping potential. For higher T or deeper BEC 
regime where phase separation is prohibited, one may use an 
rf pulse with a narrow cross section (as has been reported ex- 
perimentally t23il ) without population imbalance to achieve 
similar results. In this way, the narrow rf beam mimics the 
effect of artificial phase separation by picking up signals only 
from the central part of the trap [24] . 
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